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1. Introduction 

The search for new phases of quantum matter is as central to modern condensed matter 
physics as the historical search for new elementary particles. In recent years, this search 



H ■ has focused increasingly on frustrated magnet systems in which conventional descriptions of 

magnetic order break down entirely. 

ZnCr204 is a prototypical example of such a frustrated magnet. Despite strong (J ~ 52 K) 
antiferromagnetic exchange interactions between neighbouring Cr^+ ions, ZnCr204 fails to or- 
der magnetically down to a temperature of T/v = 12.5 K, just few percent of the Tmf ~ 400 K 
anticipated from mean field theory.^' The unconventional nature of magnetism in ZnCr204 
raises the possibility that it might also exhibit new forms of order in magnetic field. Indeed, 
it has long been known that the interplay between magnetic field and frustrated exchange 
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interactions can lead to both magnetization plateau — magnetic "solids" which break the 
translational symmetries of the lattice, and spin-flopped phases — magnetic "superfluids" 
with well-defined magnetic order in the plane perpendicular to the magnetic field. One partic- 
ularly intriguing proposal is that magnetic field might stabilize phases which break both sets 
of symmetries simultaneously. Such a phase would be a magnetic analogue of the "supersolid" 
long sought in ^He.^^ 

For a long time, testing these predictions for materials like ZnCr204 was rendered im- 
possible by the need for multiple extremes — magnetic fields B in excess of 100 T and 
temperatures of a few Kelvin. Now, a new generation of pulsed high-field magnets and ad- 
vances in instrumentation make it possible to explore this physics for the first time. Here we 
report the first determination of the magnetic phase diagram of ZnCr204, through Faraday 
rotation, for magnetic fields of up to 190 T, and temperatures down to 4.2 K. Wc find evidence 
for three distinct magnetic phase transitions, the 1/2 plateau and a 3:1 canted phase, a 2:1:1 
canted phase, all of which are well-described by a simple model based on spin-lattice coupling. 
The 2:1:1 canted phase is found for the first time in ZnCr204 as a series of chromium spinel 
oxides. This phase has been predicted by a theory of Penc et al. to realize in a small limit of 
the spin- lattice coupling. The new 2:1:1 phase is regarded as the reentrant supersolid phase 
after the 3:1 phase appeared in this system, according to an exact magnetic analogy proposed 
by Matsuda and Tstuneto,''^ and Liu and Fisher. '^^ The finite temperature evolution of these 
phases is studied, and compared with Monte Carlo simulations of an effective spin model. 
These results provide a beautiful illustration of the phenomenon of "order by distortion" 
the subtle interplay between spin and lattice degrees of freedom in a frustrated magnet. 

2. Theoretical Background 

The spinel oxide, ZnCr204 is an electrical insulator which bears all the hallmarks of a 
frustrated magnet. The metallic Cr ions are in a 3-|- ionization state, with a magnetic moment 
of 3 fiB, and occupy the sites of a highly-frustrated pyrochlore lattice, built of corner-sharing 
tetrahedra (cf. Fig. 1). At high temperatures, the magnetic susceptibility of ZnCr204 shows 
Curie- Weiss law behaviour, with a Curie temperature of 9cw ~ 390 K,"^'^^ suggesting that 
neighbouring Cr^"^ ions have antiferromagnetic exchange interactions J « 52 K. However in 
the absence of magnetic field, ZnCr204 does not order magnetically unless cooled to T^r = 12.5 
K, just a few percent of Ocw-^'"^^ The onset of this magnetic order is accompanied by a 
distortion of the lattice from cubic to tetragonal.^) 

The figure of merit Tn/9cw ~ 0.03, suggests that ZnCr204 is one of best known approxi- 
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Fig. 1. (Color online) (a) Pyrochlore lattice of magnetic Cr'^+ ions within the spinel oxide ZnCr204, 
shown in the cubic, high-temperature, unit cell. Cr"^"*" ions with magnetic moment 3 /is are situated 
at the corners of each tetrahedron, (b) Simple model of spin-lattice coupling in Cr spinels : anti- 
ferromagnetic bonds with Si • Sj < contract, while ferromagnetic bonds with Si • Sj > extend, 
in order to maximize the energy gain from exchange interactions. This systematic distortion of 
individual bonds lifts the frustration inherent in the pyrochlore lattice, stabilizing unconventional 
magnetic ground states. 

mations to an ideal (classical) frustrated antiferromagnet with a pyrochlore lattice possessing 
a three dimensional antiferomagnetic nearest neighbor interaction, for which = 0, and 
hence Tn/Oqw = 0.^'^^ It is therefore a very attractive system to study in magnetic field. 
Such high-field experiments have already been performed for the sister compounds HgCr204^^ 
and CdCr204.^^ These show a dramatic half-magnetization plateau with m = 3/2 fis-'^'^^^^^ 
Again, the transition into the half-magnetization plateau is accompanied by a structural 
distortion of the lattice. 

It is clear that magnetic order in the ACr204 spinels is intimately linked to lattice struc- 
ture. In fact the main features of the magnetic phase diagrams of IIgCr204 and CdCr204 are 
well-described by a simple theory due to Penc et al, which describes the spin-lattice coupling 
within the tetrahedra making up the pyrochlore lattice, ^^"^^^ see also ref. 16. Since exchange 
interactions originate in the overlap of electronic orbitals, magnetic energy can be gained 
by shortening antiferromagnetic bonds and lengthening ferromagnetic ones (cf. Fig. 1). This 
effect, known as magnetostriction, plays a secondary role in conventional magnets. However 
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it is essential to understanding Cr spinels, where it lifts the degeneracy associated with the 
pyrochlore lattice, and has been dubbed "order by distortion".^) 

In its simplest form, magnetostriction can be accounted for by an effective, classical, spin 
model 

-H = jY,[S^-Sj-b{S^-Sjf]-hY,S!, (1) 

where spin-lattice coupling is parameterized by the biquadratic interaction b, (ij) counts the 
nearest-neighbour bonds of a pyrochlore lattice, and h is the external magnetic field, measured 
in natural units.^' This effective spin model can easily be generalized to take into account 
longer range exchange interactions, and has the great advantage that it is accessible to Monte 
Carlo simulation. The different ordered states found can be classified according to the Ai, 
T2 and E irreducible representations of the tetrahedral symmetry group Td.^^'^^) For large 
values of b, eq. (1) favours collinear states, including a collinear half-magnetization plateau, 
in which three spins point up, and one down, within each tetrahedron (a state with local T2 
symmetry). -^^^ 

Empirically, the strength of the spin-lattice coupling b is found to depend on the metallic 
A-site ion. Experiments on HgCr204 and CdCr204, for which spin-lattice coupling is strong, 
reveal dramatic half-magnetization plateaux and magnetic phases diagrams in strikingly good 
agreement with simulations of eq. (1), for values of 6 f« 0.1.^^^ Experiments under pressure con- 
firm that &cd < ^Hg) and suggest that spin-lattice coupling is weaker again in ZnCr204.^''''^^) 

For small values of b, eq. (1) predicts a still-richer phase diagram, dominated by coplanar 
phases which simultaneously break spin-rotation symmetry in the plane perpendicular to 
the magnetic field, and exhibit modulation in the component of the spin, aligned with it. 
These phases are a natural generalization of the magnetic supersolids proposed by Matsuda 
and Tsuneto^) and Liu and Fisher.^) These and subsequent authors investigated magnetic 
supersolid stabilized by the combination of magnetic field and easy-axis anisotropy. However 
the theory of Penc et al}^^ suggests that the combination of magnetic field and spin- lattice 
coupling provides a new route to supersolid phases, even in the absence of magnetic anisotropy. 
In particular, for small values of spin-lattice coupling, it predicts an unusual 2:1:1 canted phase 
bordering on the half-magnetization plateau, which is a magnetic supersolid. This state is not 
present for the larger values of spin-lattice coupling realized in HgCr204 and CdCr204, and 
has yet to be observed in experiment. Could this phase be found in ZnCr204? 
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Fig. 2. (Color online) Schematic plan of experiments to measure the magnetization of ZnCr204 in 
high magnetic field. Extremely high magnetic fields are generated by discharging a bank of capaci- 
tors through a single-turn coil. The sample is illuminated with light from a laser, linearly polarized 
at 45° to the vertical axis. The polarization of this light undergoes Faraday rotation through an 
angle which depends on the magnetization of the sample. Horizontally and vertically polarized 
components of the transmitted light are resolved, and used to reconstruct the magnetization of the 
sample during the magnetic field pulse. The sample is held at a constant temperature of between 
4.2 and 24 K by a liquid- He flow type cryostat (not shown). 



3. Experimental Procedures 

Measurements of ZnCr204 {T^ /Ocw ~ 0.03) in magnetic field are clearly desirable, since 
it is a better approximation to an ideal frustrated pyrochlore antiferromagnet than either 
HgCr204 {Tn/Ocw ~ 0.18)2) or CdCr204 (Tn/Ocw ~ 0.11),i°) and posses a relatively small 
spin-lattice coupling. However the scale of the antiferromagnetic interactions in ZnCr204 im- 
mediately present a problem, since the saturation field for an ideal pyrochlore antiferromagnet 
with J ~ 52 K is of order 400 T — well beyond the reach of any conventional magnet. Such 
fields are however accessible through electro-magnetic flux compression techniques developed 
recently at the Institute for Solid State Physics of the University of Tokyo, which can generate 
up to 700 T.^^) Moreover, the low value of the Neel temperature T/v = 12.5 K in the absence of 
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Fig. 3. (Color online) Profile of magnetic field (lower panel), and intensity of transmitted light (upper 
panel) as a function of time, during a 190 T magnetic field pulse generated by a single-turn coil. 
Light is resolved into vertical (s) and horizontal (p) components with intensity 1^ and Ip, and the 
angle of the rotation of the light is then given as 6*^ = -i sin^^( j^^p). Since the Faraday angle 
9f is directly proportional to the magnetization m of the sample, this leads to the magnetization 
curve shown in Fig. 4. The magnetic field curve was obtained by integrating the voltage induced in 
a separate pick-up coil. The backdrop to these plots shows the solid copper single-turn coil before 
and after the experiment. 



field suggests that experiments must be performed at temperature of a few Kelvin. To achieve 
these multiple extremes, we have used a hand-made liquid-He flow type cryostat made totally 
of a "stycast" resin, ^"^^ which fits within the single-turn copper coil of a pulsed magnet, capa- 
ble of producing magnetic fields approaching 200 T for a duration of ~ 10 /xs.^^'^^-* In these 
experiments, the inner diameter of the single-turn coil was 10 mm. A current of the order 
of 2-3 MA was injected into this coil by discharging 200 kJ fast capacitor banks charged to 
50 kV. The value of magnetic field during the pulse was measured by a calibrated pick-up coil 
located in the vicinity of the sample as shown in Fig. 2. The calibration of a pick-up coil is 
described in detail in ref. 22. The estimated error of the absolute value of the magnetic field 
found by this method is about ±3 %. The lower panel of Fig. 3 shows the waveform of a 190 T 
magnetic field pulse generated by the single-turn coil method. 

The remaining challenge is to determine how the magnetic properties of ZnCr204 change 



6/13 



J. Phys. Soc. Jpn. Full Paper 

with magnetic field, within the microsecond timeframe of the pulsed field. Under such extreme 
conditions, using a traditional electro-magnetic induction method to measure the magneti- 
zation of the sample has many associated difficulties. The greatest of these arises from the 
huge background voltage induced in the pick-up coil.^^'^^^ This problem can be avoided by 
measuring the magnetization optically, using the Faraday rotation method. The experimental 
setup is shown schematically in Fig. 2. A semiconductor laser (a coherent "Cube") of the 
wavelength 635 nm was used as a light source. The light incident on the sample was linearly 
polarized at 45° to the vertical ax;is. The transmitted light was separated into perpendicular, 
linearly polarized s and p components by a Wollaston prism, as shown in Fig. 2. Raw data for 
the intensities Ig and Ip of the transmitted light are shown in Fig. 3. The magnetization of the 
sample was calculated from these as M (x Op = ^sin^^ {yTT^- Independent measurements 
in a long-pulse magnet for fields of up to 50 T demonstrate that the induced magnetization 
of the sample can be accurately measured by this method, once appropriate allowance has 
been made for the diamagnetism of the sample and quartz substrate.^) In this earlier work, 
we have demonstrated that the Faraday rotation method correctly reproduces independent 
measurements of the sample magnetization for fields of up to 50 T, and used it to measure 
the complete magnetization curve of CdCr204 in magnetic fields of up to 140 T.^) 

A single crystal sample of ZnCr204 was grown by a vapour transport technique. The 
sample was cut parallel to the (111) crystal surface, attached on a f|uartz substrate and 
polished to about 100 /^m thickness. The diameter of the resulting sample was about 1.5 mm. 
The magnetic field was applied parallel to the [111] direction. This direction is equivalent in 
all existing domains of ZnCr204, because there are three domains in which hard axes are 
along [100], [010], or [001] directions below Tjv-^^'^^^ 

Monte Carlo simulations were performed using a local-update Metropolis algorithm to 
sample spin configurations. We typically perform 10^ MC samplings for measurements after 
10^ steps for thermalization. We have checked the convergence by comparing the results for 
different initial spin configurations. In particular, to minimize the hysteresis associated with 
first-order transitions, we used mixed initial conditions in which different parts of the system 
are assigned different ordered or disordered states. 

4. Results and Discussion 

4-1 Identification of novel magnetic supersolid phase 

In Fig. 4(a), we present results for the magnetization of ZnCr204 at a temperature of 
9 K, for fields B of up to 190 T. The magnetic phase transitions identified from this data 
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Fig. 4. (Color online) (a) The magnetization of ZnCr204 in a magnetic field of up to 190 T, as 
determined by Faraday rotation, at a temperature of 9 K (solid red line). The Faraday rotation 
9m has been corrected by a constant offset 9m — 9p ~ 59dia., to allow for the diamagnetism of 
the sample and quartz substrate. A magnetisation in absolute units [/XB/Cr'^+] is obtained by 
normalizing to an independent measurement made in a long-pulse magnet for fields of up to 50 T, 
at a temperature of 4.2 K (black dashed line). Magnetic phase transitions inferred from this data 
are labelled with Greek characters (a, /3, and 7). (b) Theoretical prediction for the magnetization 
of a Cr spinel with antiferrimagnetic nearest neighbour exchange J, and weak spin-lattice coupling 
b = 0.02, in applied magnetic field, following ref. 12. Spin-lattice coupling relieves the frustration 
associated with the pyrochlore lattice, stabilizing new magnetic phases in applied magnetic field. 
The spin configurations within a single tetrahedron are shown for each phase. 



are labelled a, /3 and 7. For B < = 120 T the magnetization is linear in field, as would 
be expected for an antiferromagnet with canted spins. In the theory of Penc et al, shown in 
Fig. 4(b) for a spin-lattice coupling b = 0.02, this phase is a 2:2 canted state in which spins 
are canted in pairs within each tetrahedron (E symmetry state). 

At Ba = 120 T the system undergoes an abrupt, first-order transition into a state with 
higher magnetization and finite, field-independent magnetic susceptibihty, labelled "Phase X" 
in Fig. 4(a). We identify this phase with the previously unobserved 2:1:1 canted phase pre- 
dicted by Penc et al. for small values of spin-lattice coupling — cf. Fig. 4(b). This unusual 
spin configuration is stabilized by the competition between spin and lattice degree of freedom, 
and mixes the T2 and E symmetries of a single tetrahedron. Moreover, since the modulation 
of the component of the canted spins breaks the space group symmetries of the lattice 
("solid" order), while their and S'^ components break spin rotation symmetry ("super- 
fluid" order), this phase is a magnetic supersolid. This novel supersolid phase persists up to 
Bjs = 135 T, when the systems undergoes a continuous phase transition in a state with vanish- 
ing magnetic susceptibility and exactly half the saturation magnetization. This is the collinear 
half-magnetization plateau with a local T2 symmetry, previously seen in CdCr204^' ^'^^ and 
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HgCr204.2) At = 158 T the syst em undergoes a second continuous phase transition into 
another phase with finite magnetic susceptibihty. We identify this with the 3:1 canted T2 state 
predicted by Pcnc et al, which is also a magnetic supersolid. 

Within the theory of Pane et al. , the ratio u = W/ of the onset field for the half- 
magnetization plateau to its width W provides a measure of the effective spin-lattice coupling 
h. For HgCr204, = 10 T and W = m (wRg = 1-7), while in CdCr204, in = 28 T and 
= 34 T (wcd = 1-2). In comparison, in ZnCr204, B^ = 135 T, and = 23 T (wzn = 0.17). 
It follows that uj-Rg > <^Cd ^ ^Zn- This is entirely consistent with observation of a 2:1:1 canted 
state and the finding from pressure experiments that the spin lattice coupling in Cr spinels 
becomes stronger as the size of the A-site ion is increased, i.e. 6zn ^ bcd < &Hg-^^'^^^ 

4-2 Evolution of phases at a finite temperature 

Measurements of the magnetization of ZnCr204 were made for a range of temperatures, 
illustrated in Fig. 5(a), and showed a strong temperature dependence. For comparison, in 
Fig. 5(b), we show the predictions of eq. (1), taken from Monte Carlo simulations with b = 
0.02. These experimental results permit to construct a temperature-field phase diagram for 
ZnCr204. This is shown in Fig. 6(a), where temperature (magnetic field) has been normalized 
to the value of J = 52 K ( J = 38.9 T).^^) 

An equivalent phase diagram can also be constructed from these Monte Carlo simulations, 
where phase transitions are identified from peaks in the relevant order parameter susceptibil- 
ities. The results for spin-lattice coupling b = 0.02 are shown in Fig. 6(b), where (arbitrary) 
four-sublattice order was enforced by the inclusion of a ferromagnetic third-neighbour inter- 
action J3 = —0.05.^^) Overall the agreement between theory and experiment is excellent, 
considering the small number of adjustable parameters used (6 and J3). Critically, the topol- 
ogy of the experimental phase diagram is correctly reproduced by the theory, with the new 
phase identified with a 2:1:1 state (a magnetic supersolid) sandwiched between the 2:2 canted 
antiferromagnet and half-magnetization plateau. The main difference between theory and 
experiment is seen in the absolute value of the transition temperatures from the param- 
agnetic phase to the antiferromagnet or half-magnetization plateau phase. However these are 
highly sensitive to longer-range interactions, which in turn will be modified by the structural 
transitions associated with each magnetic phase transition — an effect not treated in the 
theory. 

It is also hard to give precise estimations of critical fields from experiment, since small 
changes in the magneto-optical data are hard to resolve at higher temperatures. However 
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Fig. 5. (Color online) (a) Magnetization of ZnCr204 obtained from Faraday rotation in fields of up 
to 190 T, for temperatures ranging from 5.5 K to 24 K. Arrows show the location of magnetic 
phase transitions. Errors in the measurement of temperature and magnetic field are estimated to 
be ±1 K and ±3 %, respectively, (b) Magnetization of bilinear-biquadratic Heisenberg model with 
b = 0.02, as a function of magnetic field h, as determined by Monte Carlo simulation. This model 
captures the essential features of spin-lattice coupling in Cr spinels. Simulations are performed for 
temperatures ranging from 0.04J to 0.32J, where the nearest neighbour exchange J « 52 K in 
ZnCr204. Arrows show the location of magnetic phase transitions. 

we observed a strong temperature dependence of critical fields, as shown in Fig. 6(a). The 
width of the half-magnetization plateau state widens as temperature increases, as a result of 
stabilization of collinear states by thermal fluctuations. This is an example of the entropy- 
driven "order by disorder" effect, ^^'^^^ and it becomes more evident at the limit of small 
spin-lattice coupling. 

5. Conclusions 

We have used Faraday rotation to establish the magnetic phase diagram of ZnCr204 in 
pulsed fields of up to 190 T, for temperatures down to 4.2 K. Four distinct magnetically ordered 
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Fig. 6. (Color online) (a) Experimental magnetic phase diagram of ZnCr204 as found by Faraday 
rotation, (b) Theoretical prediction for magnetic phase diagram of ZnCr2 04, as determined by 
Monte Carlo simulation of a bilinear-biquadratic Heisenberg model capturing the essential features 
of spin-lattice coupling. The parameter for this model are the same as those used to calculate the 
magnetization in Fig. 5(b). 

phases can be distinguished at low temperatures; these include a half-magnetization plateau, 
and three phases with canted spins and a finite magnetic susceptibility. The key features of 
these results are all well-described by a simple model based on spin-lattice coupling. On the 
basis of this analysis, we find that spin-lattice coupling in ZnCr204 stabilizes two canted 
phases bordering on the half-magnetization plateau, both of which are magnetic supersolids 
in the sense of Matsuda and Tsutneto or Liu and Fisher. One of these phases is an unusual 
2:1:1 canted state, predicted to occur in the limit of small spin-lattice coupling, but never 
previously observed. We conclude that the combination of magnetic field and spin-lattice 
coupling can stabilize novel magnetic supersolid phases in frustrated magnets, even in the 
absence of easy-axis anisotropy. 
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